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ABSTRACT 


A numerical study of the response of aircraft wings to atmospheric 
gusts and to nuclear explosions when flying at subsonic speeds is presented. 
The method is based upon unsteady quasi-vortex-lattice method, unsteady 
suction analogy and Pad^ approximant. The calculated results, showing 
vortex lag effect, yield reasonable agreement with experimental data for 
incremental lift on wings in gust penetration and due to nuclear blast 
waves . 
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NOMENCLATURE 


At 

m 

Bi 


c* 

C{k) 

Cc(k) 

C] 

Clo 

F 


= coefficients of Pad^ approximant 
= wing aspect ratio 

= coefficients of partial fraction of Padd approximant 
= local span at the section containing the transducer as 
shown in Fig. 17(b) 

= reference chord 
= semi chord length of the airfoil 
= Theodorsen function 
= generalized Theodorsen function 
= two-dimensional oscillatory lift coefficient 
= two-dimensional steady state lift coefficient 
= oscillatory pitching moment due to sinusoidal gusts 
= Fourier transform 


F(k) 

h(s) 

Jo(k), Ji(k) 

k 

k' 


L 

L 

L'(k) 

L"(k) 

M 


r 

s 
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= C(k)[Jo(k)-iJi(k)] + iJi(k) 

= impulse lift function due to nuclear blast 
= Bessel functions of the first kind 
= oic/U , reduced frequency 
= k/0.61, effective reduced frequency 
= Laplace transform 
= oscillatory lift distribution 
= in-phase component of oscillatory lift coefficient 
= out-of-phase component of oscillatory lift coefficient 
= freestream Mach number 
- ik, Laplace transform variable 
= Ut/c* or Ut/c, non-dimensional distance parameter 
= reference wing area 
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T(k) 

U 

V(k) + iW(k) 
Wa 


Wa* 
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wg* 
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= time coordinate 

= L' (k) + iL"(k) , oscillatory lift coefficient due to 
sinusoidal gusts 
= frees tream velocity 
= l-T(k), [see Eq. (30)1 

= vertical velocity component of the vortex sheet simulat- 
ing the airfoil in gust field 
= amplitude of 

= Wa’^/U 

= vertical velocity component of gust 
= amplitude of vertical velocity component of gust 

= Wg*/U 

= vertical velocity component of a sharp-edged gust 
= x/c* or x/c, non-dimensional x-coordinate 
= x-coordinate of wing leading edge 
= reference point for gust phase 

= distance penetrated into a step gust, [see Fig. 17(b)] 
= Cartesian coordinate system [see Fig. 1] 

= angle of attack 
= 


3i 


5(t) 

e 


p 

^(s) 


= root of the polynomial in the denomenator of Pad^ 

approximant, which is always real 
= delta function 

= blast intercept angle (see Fig. 20) 

= frees tream density 

= phase difference (in degrees) of oscillatory lift 
= indicial lift due to step gusts 
= radian frequency associated with the gust wavelength 
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1. INTRODUCTION 


Estimation of the response of an aircraft due to atmospheric gusts 
has been the subject of numerous investigations from the viewpoint of 
producing useful data on the induced aerodynamic forces for the design of 
active control systems for gust alleviation. 

In theorectical analyses, the change of lift and moment on a wing 
passing through a sharp-edged gust was first calculated for incompressible 
two-dimensional flow by Kdrmcin and Sears with simple mathematical formulae 
(Ref. 1). Miles (Ref. 2) extended the calculations to a travelling gust 
field, i.e., sharp-edged gust moving either downstream or upstream rela- 
tive to the airfoil. Drischler and Diederich (Ref. 3) presented results 
for a wide range of wings in both incompressible and compressible flows. 
Meanwhile, the response of an ai; foil entering a hannonic gust field was 
first introduced by Sears (Ref. 4). Murrow, et al . (Ref. 5) provided many 
numerical results of lift and moment for finite wings moving through a 
harmonic gust. Giesing, et al . (Ref. 6) also furnished some good sug- 
gestions in computing the oscillatory lift and moment. One notable method, 
called the Doublet Lattice Method (DLM) which was originally developed by 
Albano and Rodden (Ref. 7), was later improved to become a very useful 
tool in unsteady aerodynamics (Refs. 8, 9). 

For the general harmonic analysis, the atmospheric gust was con- 
sidered as a random set of discrete gusts. Response had been predicted 
most commonly with the assumption that the vertical component of gust 
varied along the flight path, but did not vary along the span. This 
assumption was adopted by most researchers because the associated compu- 
tations were less extensive than those for the more general cases of 
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random gusts. It may not be sufficiently accurate for very large aircraft, 
but It should provide useful data for most configurations. 

Besides all these numerical calculations, not much experimental work 
appeared to have been done or to be available. Roberts and Hunt (Refs. 10, 
11) made a series of measurements of transient pressures on a narrow delta 
wing of M"1.2 due to vertical gusts, and Patel presented some experi- 
mental results for a couple of delta wings (Ref. 12) and other types of 
wings (Ref. 13) In harmonic gust fields. 

In current aerodynamic research, the vortex flow phenomenon Is 
drawing much attention because It offers significant contributions to 
aerodynamic characteristics on low-aspect- ratio wings with sharp or thin 
edges. For these types of configurations, the pressure distributions due 
to the leading-edge vortex separation are drastically different from those 
given by the conventional linear theory. The complex flowfleld also makes 
it more difficult to predict aerodynamic forces accurately. Most of the 
aforementioned unsteady aerodynamic methods are based on the linear theory, 
and can not predict the leading-edge vortex effect. Although Atta, et al . , 
(Ref. 14) developed an unsteady lifting-surface method with vortex flow 
by using the unsteady vortex-lattice method, no application to various 
gust problems has been presented and the computing cost is expected to be 
very expensive. 

In this report, an unsteady lifting-surface method from reference 15 
Will be used to calculate the lift and pitching moment due to sinusoidal 
gusts for several wing planforms. The calculated results are compared with 
other theories for the attached flow and with experimental data for vortex 
flow. In calculating the response of a wing to a sharp-edged gust, Garrick 
(Ref. 16) developed a relation between the oscillatory forces due to flight 
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through continuous harmonic gusts and the indicia! forces due to sharp- 
edged gusts. Instead of Fourier transform taken by most other theories 
in handling this reciprocal relationship, the present method will use Padd 
approximant to represent the harmonic response and Laplace transform will 
be used to calculat<^ the indicial functions. The problem formulation and 
computed results are presented in the following chapters. 

For many years, military personnel have been continously interested 
in the prediction of the response of an aircraft resulting from a nuclear 
blast wave. Kaman AviDyne (Ref. 17) did a series of experiments to mea- 
sure the blast pressures on a rigid highly sweptback wing at high subsonic 
speeds. McGrew, et al . (Refs. 18, 19) recently used the DLM to develop a 
nuclear blast response computer program for wing-body configurations. Yet 
the method is valid only for the attached flow. The present method, based 
on the unsteady suction analogy, will be used to demonstrate the capabili- 
ty of predicting the nuclear blast response involving the vortex flow. 

Some results- in simulating experimental data in reference 17 will be shown. 


2. MATHEMATICAL FORMULATION 


2.1 Two-Dimensional Gust Penetration 

Consider that a thin airfoil moving with a velocity U, enters a 
region of atmospheric gust with velocity distribution Wg normal to the 
direction of motion (see Fig. 1). The boundary condition requires that 
the total vertical velocity due to the gust and the vortex sheet simulat- 
ing the airfoil must vanish; 

Wg + Wq - 0; for z = 0, -c* < x < (1) 

where Wg is the vertical velocity component of the vortex sheet, x,z are 
the coordinate systems attached to the airfoil and c* is the semi chord 
le^‘"^th of the airfoil. 

To solve a simple harmonic gust problem, the following expression 
is used to specify the sinusoidal gust: 

Wg(x,t) = , (2) 

where Wg* is the amplitude of Wg and oj is the radian frequency associated 
With the gust wavelength. Substituting into Eq. (1) and canceling out 
the time factor Eq, (2) leads to 

Wa*(x) = - (3) 

By way of introducing the reduced frequency, k = wc*/U, Eq. (3) then be- 
comes 



g-ik(x/c*)^ 


(4) 
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It ir convenient to divide both sides by U and use the expression x-x/c* 
as the non-dimensional x-coordinate, 


U ” " U 


■j.e. , Wg = " Wg e " 

where Wg = Wg* / U and Wg * Wg* / U . 


(5) 


Using Eq. (5), the exact lift distribution for incompressible flow 
can be shown to be (Ref. 20), 

L = 2TrpUc*Wg{C(k)[Jo(l<)-iJi(k)] + iJi(k)}e^^^, (6) 

where C(k) is the Theodorsen's function, Jn(k) and Ji(k) are Bessel func- 
tions of the first kind. For compressible flow, C(k) will be replaced by 
generalized Theodorsen's function** CQ(k) and all other terms in Eq. (6) 
remain the same. 

Furthermore, the lift caused by an arbitrary Wg can be calculated 
from Eq, (6). For a sharp-edged gust striking the leading edge of the 
airfoil at t = 0, the boundary condition is 



X > Ut - c*. 

(7) 

X < Ut “ c*. 


**The generalized Theodorsen's function was generated by Mr. Chung-Hao Hsu 
in his earlier work for calculating lift on wings due to step change in 
angle of attack. 


5 



ORIGlN;-,i, : i 

OF POOR QUALnV 

Then, Wg can be represented by the Fourier integral, 


vjq = — f f(w) e dw, (8) 

^ 27T J -00 


which can be inverted into 


r ~ 

f(w) - Wg e dt. (9) 

-00 

Since Wq = 0 for t , Eq. (9) can be shown to be, by following 

y U 

the Fourier transform of a constant (Ref. 20, p.287), 


f(03) 



Woe“i‘'^(x+c*)/U 

ioj 


( 10 ) 


Substituting Eq. (10) back into Eq. (8), the boundary condition for the 
sharp-edged gust becomes 


Wn = 



i'co(t - 


u 


doj 


1(0 



gik(s - X - 1) 
ik 


dk, 


( 11 ) 


v^here s = Ut/c* is the non-dimensional distance parameter. The airfoil 
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lift due to the harmonic gust (Eq. (5)) is given by Eq. (6). Therefore, 
for the step gust based on Eq. (11), the lift can be calculated as 


L - pUc*Wo 


{C(k)[Jo(k)-iJi(k)] + iJi(k)> 


ik 


dk. (12) 


From Eq. (12), the non-dimensional lift development due to a step 
gust, i|^(s), is given by 


with 


L = 27 TpUc*Wo»|^(s) , 




2it J « 


{C(k)[Jo'(k)-iJi(k)] + iJi(k)} 
ik 


(13) 


dk. (14) 


Let F(k) = C(k) [Jo(k) - iJi(k)] + iJi(k), then 


\p(s) = 


_ 1 


F(k)e 


ik(s-l) 


2tt 


ik 


■ dk 


JL r “ 

2tt / 




ik 


(15) 


As k is always greater than or equal to zero, Eq. (15) can easily be 
inverted to 


F(k) e 


-ik 


ik 


\jj(s) e ds, 


or 


F(k) e"’''^ = ik j i|j(s) e"''’^^ 


ds. 


(16) 
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Let r= ik. Then. 
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F(k)e“'' = r 



= r L { i(/(s) }, 


(17) 


i|^(s) = L"' (18) 

where L { i|;(s) } is the Laplace transform of i|;(s). 

Hence, the indicia! lift function can be obtained from the inverse 
Laplace transform involving the amplitude of lift distribution due to a 
sinusoidal gust. 


2.2 Three-Dimensional Gust Penetration 

Now consider a rigid thin wing travelling at speed U through an 
infinite array of harmonic gusts with vertical velocities Wg, uniformly 
across the wing span. The boundary condition is similar to that for a 
two-dimensional sinusoidal gust: 

Wa = - (19) 

where k again is the reduced frequency with a reference chord length c, 
and is a reference point for the gust phase. 

Following reference 21, the oscillatory lift force due to the har- 
monic gust (Eq. (19)) is 

L(t) = pUwgS{L'(k) + iL"(k) (20) 
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OF r ... 

where S is the reference area of the wing, and L'{k), L"(k) are the in- 
phase and out-of-phase components of the dimensionless lift force. 

Next consider a wing entering a sharp-edged gust. Again the 
boundary condition is very similar to Eq. (7), 


f 0. 





Wo, 


X > Ut + x^g, 

( 21 ) 

X < Ut + x-|g, 


where X]g is the x-coordinate of the leading-edge of the wing. From Eq. 
(20), one can determine the indicial lift function representing dimen- 
sionless lift development due to a step gust, »|j(s), 


L(s) = pUwoSifj (s), 


( 22 ) 


with 


i{i(s) 


J_ f ” ( l-'0<) ■■■ il-"(k) 

2tt / ik 

“00 


(23) 


Let 

T(k) = L'(k) + iL" (k), 


(24) 


’l'(s) 



Me ^'^=dk. 
ik 


(25) 


With the same procedures from Eqs. (15) - (18), i|i(s) can be obtained from 
the inverse Laplace transform involving the dimensionless lift distribu- 
tion due to a harmonic gust, 

<p(s) = {M. }. (26) 

r 
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To be able to calculate the inverse Laplace transform indicated in 
Eqs. (18) and (26), it is convenient to express T(k) and F(k)e“*^ as closed- 
form functions of r. In any lifting-surface computation, T(k) or F(k) are 
calculated only at a finite number of k's. These values can then be in- 
terpolated by Fade approximant as suggested by Vepa (Ref. 22). Following 
Vepa, an [N,N ] sequence of Fade approximant to approximate T(k) in three- 
dimensional case and F(k)e“'^ in two-dimensional gust can be written as, 

T(k) = 1 - [ N,N ), F(k)e"*' = 1 - [ N,N 1, (27) 

[ N.N 1 = fl.rN H- A.rN-1 + + A^.ir _ (23) 

r^l + A\|r^”^ + ’••••“•• + 

where A-{ are the coefficients of Fade approximant. 

For N = 3 and use T(k) as an example, 

T(k) = 1 - Azr _ (29a) 

r®+ Aar^+ + As 

and for N = 2, 


T(k) = 1 . (29b) 

r^+ Aar + As 

Let 

1 - T(k) = V(k) + iW(k). (30) 

From Eq. (30), it follows that for N = 3, 
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V(k) + iW(k)= . (31) 

(ik)H A3(ik)^+ A4(ik)+As , 

Eq. (31) can be expanded and separated into real and imaginary 
parts respectively to give, 
real part: 

Aik^- AaVk- A4Wk + A 5 V + Wk^= 0, (32a) 

imaginary part: 

Aok^- Aak- A3Wk + A4Wk + AsW- Vk^= 0. (32b) 

Normally, more values of k's are chosen to calculate T(k) than the 
number of unknowns Ai in Eq. (32a) and (32b). Therefore, a least square 
technique to be described next is used to determine A^-'s. 

2.4 Least Square Technique 

The least square principle is based on the requirements that Aj's 
are determined by minimizing the sum of squares of errors: 

M 

Sum= ^ 

1 1 

where M is the number of k's of which T(k) is calculated by any existing 
lifting-surface theory. 

At a minimum, all the partial derivatives with respect to A-j's, 
such as 8Sum/3Ao, 3Sum/9Ai, , 8Sum/3A-| , must vanish. These condi- 

tions result in i+1 equations for i+1 unknown coefficients A-j's. Thus, 
Aj's can be solved from Eq. (33) with M different lift values at corre- 
sponding k's. 


(L.H.S. of Eq. (32a))^+ ^ (L.H.S. of Eq. (32b))^. (33) 
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2.5 Indicial Lift Function POOR 

After all coefficients of the Padd approximant are determined, Eq. 
(29a) can be rewritten as, 


T{k) = 1 - r 


Aor^+ Air + A2 

MPIWIIIII. • 

r^+ Asr^+ Aur + As 


By partial fraction method, Eq. (34) leads to 


T(k) _ J_ f Bi 

r r . 4 ^ r - tij ’ 


(34) 


(35) 


where 3i is the ith root of the polynomial in the denomenator of Eq. (34) 
and Bt is the corresponding coefficient of partial fraction in Eq. (35). 

Based on Eq. (26), the indicial lift function ip(s) can be obtained 
by applying the inverse Laplace transform to Eq. (35), 


ilj(s} 




(36) 


2.6 Nuclear Blast Response 

Calculation of the lift development of a thin wing encountering a 
nuclear blast wave will follow the same way as in computing the sharp- 
edged gust response. No major change has to be made except the boundary 
condition. The shock wave induced by the nuclear blast is assumed travel- 
ing at sonic speed. Thus, an impulse function 6(s) is used instead of the 
unit step function in the gust response condition (Ref. 24), 

Wg = Wo5(t 
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where 0 is the magnitude of the vector sum of the shock wave velocity and 
the frees tream velocity. Eq. (9) now becomes, 


/ OO 

Wo(S(t e “^'^^dt 

-OO U 


"Woe 


-io)( 


X 



(38) 


Comparing Eq. (38) with Eq. (10), the main difference is the factor "iw" 
in the denomenator of Eq. (10). This also follows from the fact that the 
impulse is equal to the time derivative of the step input. Hence, 


F( h(t) } = iwF{ A(t) } , 


(39) 


where h(t) is the unit impulse, A(t) is the step function and F{ } is 
the Fourier transform. 

Rewriting Eq. (25) for the impulse response, it follows that. 


h(s) = 2 ^/ T(k)e’*'=dk, (40) 

or, 

h(s) = { T(k) } . (41) 

The Pad^ approximant used for the indicial response analysis re- 
mains applicable, except the inverse transform calculation has changed. 

For example, N = 2 Fade approximant can be written as, 

T(k) = 1 - h'tLkL- 
r^+ Aar + As 
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(Ai - AoA 2 )r - A 0 A 3 
+ A2T + A 3 


9 

V 


= 1 - Ao 


2 

Z 

i==l 



From Eq. (41), the impulse lift function h(s) can obtained, 


(42) 


h(s) = { T(k) } 


= (1 - A )6(s) - 


2 

Z B-j e 
i=l 



(43) 
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3. NUMERICAL RESULTS AND DISCUSSIONS 


3.1. Gust Penetration 

3.1,1. Sinusoidal gust problem 

The lift development due to a harmonic gust is calculated by the 
computer program based upon the unsteady quasi-vortex-lattice method de- 
veloped by Lan (Ref. 23). 

For a thin airfoil, 30 vortex elements are used in the computation. 

The steady state two-dimensional lift values are simply calculated by the 

equation Ci = 27 t/ 3, where 3 =/l-M^, M is the freestream Mach number. The 
'0 

computed results are compared with Sears' for incompressible flow and with 
Graham's for several different Mach numbers at reduced frequency 2. Both 
comparisons, showing good agreement, are tabulated in Tables I and II . 

The three-dimensional unsteady aerodynamic program of reference 15 
is then revised to account for the gust response. The present attached 
flow results of lift and moment for a delta wing of 75° sweep at M = 0.4 
are compared in Figs. 2 and 3 with those calculated by a kernel function 
method (Ref. 5). 

Results of calculation will also be compared with experimental data. 
In references 12 and 13, a gust tunnel which could generate a sinusoidal 
vertical gust was used to measure the oscillatory lift and moment on two 
delta wings of AR=1 and 2, and several other commonly used wing planforms. 
The tests were performed for all wings at two mean freestream velocities 
of 12.43 and 20.00 m/s. However, the oscillatory gust waves convected 
downstream with a velocity of 0.61 of the mean freestream velocity. This 
would indeed influence the gust wavelength, i.e., the frequency parameter. 
Thus, in the present calculations based on gust moving with the freestream 


15 


velocity, an effective frequency (k' »k/0.61), as suggested by Patel, will g 

i 

be used 1n the following comparisons. 

In Fig, 4, test data for a rectangular wing of PR ~ 6 (Ref. 13) are 
compared with two sets of theoretical results. It is seen that the pre- 
sent theory agrees well with Graham's (Ref. 24} in the predicted in-phase 
component of the oscillatory lift. phase lag is underpredicted and 
both theories overpredict the lift amplitude. 

To demonstrate the vortex flow effect, a delta wing of PR=1 is used 
to compare with Patel's data (Ref. 12) at a=0 and 12 degrees in Fig. 5. 

The present theory predicts well the lift amplitude at both angles of at- 
tack. At a=12°, the phase lag is also adequately predicted by the pre- 
sent theory. On the other hand, the phase lag at a = 0° is not accurately 
predicted. In seeking the reasons for these deviations, several points 
should be noted: 

(1) At a=0 degree, the attached flow prevails. It is of interest 
to compare the present results with the doublet-lattice method (Refs. 8,9) 
for this delta wing of PR=1 . In Fig. 6, agreement between the present 
results and doublet-lattice method's is excellent in lift amplitude and 
phase angles. 

(2) As depicted in reference 12, force measurements were made re- 
lative to the undisturbed freestream gust at the root 2/3rd chord point. 

Another point at the gust tunnel exit was also used as a reference. It 
is not known whether the conditions at the exit point were disturbed once 
the test model was placed in the tunnel. 

(3) Also, Patel indicated in reference 12 from test data that the 
incremental lift due to vortex lift contribution was important in magni- 
tude only with no measurable contribution to phase angles. However, the 
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present theory shows that this is approximately true only with respect to 
some reference points Ci ,e. , X|^ at Eq. (19)1. This is illustrated in Fig. 
7 with the root midchord point and wing apex as reference points, There- 
fore, the results by the present theory very much depend on the precise 
location of the reference about which the phase angle is calculated, 

At any rate, reasons for the discrepancy in the predicted phase 
angles at a = 0 degree for the delta wing of PR ^ I are still unknown at 
the present time. 

3.1.2. Fade approximant 

Two different sequences of Pad^ approximant are constructed here 
to fit various freestream conditions. For airfoils in compressible flows, 
Pad6 approximant with N=3 in Eq. (27) (called Pad^ A6) is used. On the 
other hand, Padd approximant with N=2 (called Pade A3) is used for air- 
foils in incompressible flow and also for three-dimensional conditions. 
These choices are made through numerical correlations. The corresponding 
matrices for deciding the unknown coefficients A-,* in Eq. (28) are present- 
ed in Appendix A. 

The good agreement between the calculated results from the unsteady 
QVLM and by Pade approximant for a thin airfoil in the harmonic gust is 
shown in Fig. 8 for incompressible flow, Fig. 9 for Mach 0,5 and Fig. 10 
for Mach 0.7. The oscillatory lift for a delta wing of = 1.2 to be used 
for calculating step gust response is shown in Figs. 11 and 12 for upward 
and downward gust respectively. 

Most other theories use Fourier transform to handle the reciprocal 
relationship between the oscillatory and indicial lift forces. So they 
need harmonic lift values at a large number of different frequencies. 
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But the number of k required by the pv'esent method 1s only 7 (Including 
k = 0), and the Indiclal lift results, which will be seen In the next sec- 
tion, still show good accuracy. 

3.1.3. Indiclal lift function 

The Indiclal lift results performed by using the inverse Laplace 
transform described in the preceding chapter will be presented. First, 
for a thin airfoil passing through step gusts In Incompressible and com- 
pressible flows, Fig. 13 shows plots for Kiissner function, and Figs. 14, 

15 exhibit the Indiclal lift at Mach 0.5 and 0.7. Good results have 
been expected because of the accurate approximant shown in last section. 
The exact solutions are calculated through the data in Table 6-2 of refer- 
ence 20, From Fig. 16, it is apparent that the compressibility effect 
decreases the rate of lift build-up in two-dimensional flow. 

Second, in thin wing gust penetration, for lack of experimental 
force data to make a direct comparison, some pressure data (Ref. 11) will 
be employed to compare the trend produced by calculated total lift. 

The configuration of the delta wing used in reference 11 is shown 
in Fig. 17(a). The model wing was carried along a straight railway track 
on a rocket-propelled sledge, through the efflux from an open- jet wind 
tunnel blowing across the track. The velocity of the sledge was 180 ft/ 
sec,, the tunnel efflux velocity was 47 ft/sec., and the model was at zero 
angle of attack. For measuremr'^ts of transient pressure, four transducers 
were positioned at locations being 0.3, 0.4, 0.5, 0.6 root chords aft of 
the apex, and along a line at 75% semispan [see Fig. 17(a)]. In accord- 
ance with figures in reference 11, the indicial lift is plotted againsc 
local distance parameters x^/b^ in comparing with pressure values through 
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transducers A, B and C, where and bp are defined in Fig. 17(b). 

Fig. 18 shows the similar trend between measured transient pressure 
and the calculated indfcial lift for both upward and downward sharp-edged 
gust. It is seen that the development of a gust-induced gain of lift is 
very gradual. On the contrary, the gust- induced loss of lift occurs 
relatively instantly. In the calculations, the vortex lag described in 
reference 15 is assumed to be present if the lift is increasing and there 
is no vortex lag if the lift is decreasing. Ihis assumption appears to be 
reasonably accurate for the leading-edge vortex flow. This can also be 
seen from Fig. 19 which illustrates the comparisons among the vortex flow 
and potential results. There is some significant difference for the up- 
ward gust while the trend is quite close in the downward gust condition. 

3,2. Nuclear Blast Response 

The unsteady aerodynamic program of reference 15 is again used to 
calculate the nuclear blast responses of aircraft flying at high subsonic 
speeds. 

Reference 17 is the only obtainable test data which can be used to 
check the leading-edge vortex separation effect, predicted by the present 
method, on aircraft nuclear blast response. The sideview of a thin wing 
intercepted by nuclear blast waves is shown in Fig. 20. Fig. 21 shows the 
test model which consists of a swept wing of 67° leading-edge sweptback 
angle with a nose and partial fuselage section. In the test, the model 
was mounted on a high speed dual rail rocket sled at an initial angle of 
attack of 3.2° . The sled, travelling at Mach 0.76, was intercepted pro- 
gressively by blast waves from sequential detonation of charges of TNT 
with the blast intercept angle 0 = 20° (see Fig. 20). Twenty pairs of 
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pressure transducers were installed on the left wing half to measure the 
blast- induced pressures. 

For the purpose of correlation with the test model, a semiwing used 
by the present method is illustrated in Fig, 22, For such a configur<?,i:ion, 
the concept of augmented lift is included in the present calculation. The 
definition for the characteristic length is adopted from reference 26. With 
negative augmented vortex lift, the vortex lift effect may not be as strong 
as expected even the leading-edge sweptback angle is 67° for the test model. 

Fig. 23, reproduced from reference 17, shows the pressure variation 
measured by transducer 13 which was positioned at half semispan and along 
quarter chord line. Like DLM, the present theory underpredict the blast- 
induced incremental pressure loadings because the nonlinear vortex effect 
is not included in the calculation of pressure differential. The present 
theory, being based on unsteady suction analogy, can only demonstrate the 
vortex effect by the variation of total lift or moment. Fig. 24 shows the 
comparison of vortex flow with attached flow for impulse lift. It is evi- 
dent that the vortex lag decreases the rate of lift decay after the blast 
intercept. 

There are several factors which must be mentioned in connection with 
Figs. 23 and 24: 

(1) The blast amplitude is determined by peak material velocity 
(gust) behind the shock front and is assumed invariant throughout the 
present calculation. 

(2) In the experiment, because of different locations of transducers, 
the blast intercept time (time at the shock arrival) is different for each 
transducer. In the present calculation, the coordinate origin is set at 
the apex of the wing, so that the intercept is exactly at s=0 as 
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shown In F1g. 24. However, to correlate with the test data, a shift has 
been made in Fig. 23. 

(3) The steady state lift value is used for all time parameters 
less than zero since there is no incremental lift or pressure differential 
values for s < 0 in the present analytical results. 

For a more significant vortex lift effect, the impulse lift on a 
delta wing of 4R = 2 at Mach 0.5 intercepted by a blast wave with the same 
intercept angle as that in Fig. 20 is considered. The initial angle of 
attack is assumed to be 15“ and the peak blast-induced angle of attack 
is assumed to be an additional 15“. The attached-flow oscillatory lift 
can be well represented by a Padg approximant as indicated in Fig. 25. 
However, with the vortex lift effect included, the Fade approximant 
fails to approximate the calculated results accurately as shown in Fig. 

26. Because of this reason, the calculation of impulse lift for this 
delta wing was not successful; and hence, the results are not presented 
here. 


4. CONCLUSIONS AND RECOMMENDATIONS 


An unsteady lifting-surface computer program based on quasi -vortex 
lattice method along with leading-edge suction analogy hac been developed 
to estimate the oscillatory air forces on wings of general planforms in 
gust flow at any frequency. Pad^ approximant and Laplace transform have 
made it practical to convert the oscillatory air forces to indicia! air 
forces . 

Both the experimental data and other theoretical results are used 
to check the accuracy of the present calculations in attached flow and 
with vortex lift effect. It is shown that the present method can accu- 
rately predict the oscillatory and indicial lift on wings in different 
gust fields. Also, the phenomenon of the gust-induced gain of lift being 
very gradual and the gust-induced loss of lift occurring relat’vely ab- 
ruptly can be explained by the presence or absence of vortex lag effect. 

The present program is extended to account for the nuclear blast 
response as well. Though there is no lift data available at this time, 
the trend for the vortex lag is clearly seen from the comparison between 
the attached flow and vortex results predicted by the present method. 

The following points should be noted to improve the efficiency 
and capability of the present method. 

(1) In the present calculations, 72 vortex elements were used for 
half wingspan -• 6 in the chordwise direction and 12 in the spanwise di- 
rection. It is recommended that 40 elements (e.g,, 4 in the chordwise 
direction, 10 in the spanwise direction) could be used in lieu of 72 for 
small reduced frequencies; thus, the size of the aerodynamic influence 
coefficient matrice could be greatly reduced. 
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(2) The present method can only deal with nuclear blast locations 
being on the plane of symmetry, i.e., the X-Y plane. For more general 
blast orientations, the present method should be extended to calculate the 
aircraft wing response of the blast waves coming from any arbitrary direc- 
tion, since an asymmetric condition can always be treated as a combination 
of a symmetric and an antisymmetrical conditions. 

(3) Among the existing approximation schemes for the unsteady aero- 
dynamic loads, Karpel's approximant (Ref. 28) provides better accuracy than 
Fade's. It is suggested the present method should include more than one 
approximant to meet different gust and blast conditions. 
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Table I. Oscillatory Lift* on an Airfoil in Sinusoidal 
Gusts at M = 0. 


Present Method 


r. I ** 

Sears ' 


k 

Cl /C-, 
Ir' lo 

^li^^^lo 

^Ir'^^lo 

Cii/Ci 

0.02 

0.96289 

-0.07484 

0.963 

-0.075 

0.04 

0.92403 

-0.11449 

0.924 

-0.114 

0.06 

0.88697 

-0.13923 

0.887 

-0.139 

0.08 

0.85266 

-0.15457 

0.852 

-0.154 

0.10 

0.82126 

-0.16348 

0.821 

-0.164 

0.20 

0.70156 

-0.15964 

0.702 

-0.160 

0.40 

0.56789 

-0.08494 

0.568 

-0.085 

0.60 

0.48837 

-0.00490 

0.488 

-0.005 

1.00 

0.36865 

0.12594 

0.369 

0.126 

2.00 

0.08158 

0.26796 

0.082 

0.268 


and Cn,. are the in-phase and out-of-phase components, 
1 r n ^ 

Sears' results are copied from reference 27. 
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Table II. Oscillatory Lift on an Airfoil in Sinusoidal 
Gust at k = 2.0. 


Present Method Graham's 


M 

^Ir 

'li 

^Ir 


0.0 

0.5126 

1.6836 

0.5125 

1.6837 

0.2 

0.6925 

1.7473 

0.6955 

1.7441 

0.4 

1.1535 

1.5960 

1.1511 

1.5953 

0.5 

1.2773 

1.3139 

1.2774 

1.3138 

0.6 

1.2138 

1.1104 

1.2139 

1.1102 

0.7 

1.2192 

1.0113 

1.2152 

1.0085 

0.8 

1.1881 

0.8921 

1.1838 

0.8897 

0.9 

1.1516 

0.7986 

1.1486 

0.7961 

resul ts 

cori^espond to 

Wg = 1 in 

Eq. (4). 



Graham's results are taken from reference 23. 
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APPENDIX A Coefficient Matrices of Pad^ Approximation with Least Square 
Techniques 


1. Padd A6 (3,3) 
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OF POOR QUALITY 


f 2k®. 0 

J 

0 it 

XJ 

Ek't 0 

J 

2kfW. Zk'iU. 
J J JO 

Ekju. -Zk^W, 
0 J J J 

2k® W. Ek?U, 
J J J J 


-2k!; 

J 





'(\o" 


'' c ^ 

2k®U, 

J J 

0 





Ai 


- 2 k|w. 

J J 

-2k^ 

J 


2 k?u. 

1 J 



As 


2k*! U. 
1 1 

V 


U 0 
0 

Ek^O|?.W|) 

< 

As 

>=< 

J J 

> 

0 

"2k? U. 
J J 

0 

2k?(Uj+Wp 

0 


A4 



-2k, W. 
J J 

-Zkj(uf.w?) 

0 

j;(Uj*+wp 


As 


0 

L 




2. Padd A3 (2,2) 
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Fig. 2 
harmoni 


Generalized oscillatory lift for a 75° delta wing due to 
c gusts with Wg = 1 at M = 0.4. 





Fig. 3 Generalized oscillatory pitching moment for a 75° delta 
wing due to harmonic gusts with Wg = 1 at M = 0.4; pitching axis 
at root midchord point. 











Amplitude 


(j)/deg. 


OF Q'i^ALfuY 



Fig. 5 Incremental lift for a delta wing of /VI = 1 in sinusoidal 
gusts for a = 0 and 12 degrees with w^= 0.0314 and reference point 
at root 2/3rd chord point. 



Fig. 6 Incremental lift for a delta wing of /iR = 1 in sinusoidal 
gusts for a = 0° with Wj^= 0.0314 and reference point at root 2/3rd 
chord point. 
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Fig, 7 Variations of phase angle due to changes of angle of attack 
and reference points for a delta wing of m = 1 in sinusoidal gusts. 
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Fig. 8 Oscillatory lift for a thin airfoil due to harmonic 
gusts with vv^, = 1 at M = 0. 
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Fig. 9 Oscillatory lift for a thin airfoil due to harmonic 
gusts with w„ = 1 at M = 0 ^ 



Fig. 10 Oscillatory lift for a thin airfoil due to hannoni 
gusts with Wr, = 1 at M=0.7 
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on from a sharp-edged gust of a thin airfoil at M--0.7 
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Fig. 17(a) Gereral configuration and transducer positions of the 
test model in reference 11. 
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Fig. 17(b) Definitions of the local span and distance parameter 
penetrated into the vertical gusts. 
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Fig. 18 Continued. 
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Fig. 18 Continued. 
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Fig. 21 General configuration of the test model used for nuclear blast 
response (taken from reference 17). 
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Fig. 22 Simplified wing planform used by the present theory 
for nuclear blast response analysis. 
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Fig. 23 Pressure variation of the test wing due to nuclear blast waves with 0 = 20° for a=3.2' 
at M=0.76 (reproduced from reference 17). 
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Fig. 24 Impulse lift function of the test wing due to nuclear blast waves with 0 = 20® for 
a= 3.2® at M=0.76. 
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